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O ; ABSTRACT 

. Wc examine the possibility that hyper-accretion onto newly born, black holes occurs in 

highly intermittent, non- asymmetric fashion favorable to gravitational wave emission 
in a neutrino cooled disk. This picture of near-hole accretion is motivated by magneto- 
rotationally induced, ultra-relativistic disk dynamics in the region of the flow bounded 
from below by the marginally bound geodesic radius r-mh- For high spin values, a 
largely coherent magnetic field in this region has the dynamical implication of compact 
mass segregation at the displacement nodes of the non-axisymmetric, MRI modes. 
I When neutrino stress competes favorably for the disk dynamical structure, the matter 

, clumps may be rather dense and sufficiently long-lived to excite the Quasi-Normal 

Ringing (a.k.a. QNR) modes of the Kerr geometry upon their in- fall. We find that such 
accretion fiow may drive bar-like, quadrupole (Z,m = 2,2) modes in nearly resonant 
CO . fashion for spin parameters a > .9. The ensuing build up in strain amplitude of the 

' undamped oscillations warrants a brisk rate of energy deposition into gravitational 

r~| . waves. A detectability assessment for the LIGO interferometers through the match 

pi i' filtering technique is given by integrating the energy fiux over a one second epoch of 

resonant hyper-accretion at 1 Mq sec~^. Thus, a 15 M© Kerr black hole spinning at 
a ~ .98 (/qnr ^ 1677 Hz), and located at 27 Mpc (e.g., GRB980425), wiU deliver a 
^ characteristic strain amplitude, /ichar — 2.2_2i, large enough to be detectable by LIGO 

^ \ II. If resonant hyper-accretion were sustainable for a longer period (or at higher rates) 

possibly associated with a second broad hump in a GRB light-curve, these objects 
could be detected by LIGO I at very low redshifts. 
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1 INTRODUCTION 

It seems only fitting for ultra-relativistic, black hole accretion to be the leading contender to energize and explode the massive 
stars associated with gamma-ray bursts and hyper-nova. One must recognize, however, that such a scenario resides on a very 
exotic front: stellar-mass black holes accreting at twelve orders of magnitude above the Eddington limit! On the other hand, 
since GRB engines are effectively hidden from view, gravitational wave emission models may constitute the most effective 
probes to such enigmatic events. Indeed, there is a great deal of impetus to address such extreme accretion scenario and 
its associated, non-standard energy deposition channels. In this paper we are chiefly concerned with the possibility that the 
dynamical structure of hyper-accreting flows is highly intermittent on large scales, with the in-fall of large mass over-densities 
leading to proliflc gravitational wave emission coincident with the gamma-ray stage. 

We envisioned the accretion disk setting following the hyper-accreting black hole models of Popham, Woosley & Fryer 
(1999, hereafter PWF). For a fiducial scenario: = 3 M0, Qgg = 0.1, and M = 0.1 M0 sec~^, these authors found 
that the onset of photo-disintegration and of neutrino cooling at radii rs70 [GM/c?\ yield a mildly advective, semi-thin, 
Hs/nA, accretion disk structure with a slightly sub-Keplerian rotation profile. An improved account of neutrino transport 
(Di Matteo et al. 2002) has found that the innermost disk portion should become optically thick to neutrinos for accretion 
rates MaO.l M© sec~^. This being the case, cooling by radial advection of energy will compete with that from local neutrino 
emission and advection will overwhelm neutrino cooling for M ^1 M0 sec~^ (Di Matteo et al. 2002). Yet, such a high accretion 
rate is rather unlikely (MacFadyen et al. 2001), so we restrict our analysis to M < 1 Mq sec~^ 

The structure of this paper is a follows. In §2, we assess the relevance of physical processes occurring in the relativistic 
region of the flow: Vmb < r < Vms. Magneto-rotationally-induced, relativistic disk dynamics is addressed in §3, where we also 
discuss the role of compressibility and estimate the size of mass over-densities from a neutrino cooled disk. In §4 we build 
an idealized, analytical model for gravitational wave emission and, in §5, estimate the detectability of the GW signal for the 
LIGO interferometers with the match filtering technique. 



2 R. A. Araya-Gochez 



2 PHYSICS INSIDE THE ISCO 

The model for gravitational wave emission outlined in §4 is based on the premise that large-scale magneto-rotational effects 
will drive the disk dynamics in the ultra-relativistic region of the flow bounded from below by the marginally bound orbit 
radius. For spin rates a > .9, the growth rate of the magneto- rotational instability ctmri (a.k.a. MRI), is faster than the 
inverse dynamical time scale, n^^, throughout the annulus rmb ^ r < Tms by a factor ~ 0[2] (see Fig 1). Thus, a small delay 
to reach the "free-fall" stage is all that is needed for magneto-rotational dynamics to take place. On the other hand, rmb 
represents the lowest bound on the location of the cusp in effective potential and is akin to a relativistic generalization of the 
inner Lagrangian point, Li, for mass transfer in close binaries (Kozlowski, Jaroszyhski & Abramowicz 1978). 

The exact location of the inner boundary depends on the source of free-energy for the flow. Standard lore holds that the 
flow draws energy from the radial gradient in angular momentum and that since dissipation of angular momentum vanishes 
at rms for a Keplerian disk, this is where the inner edge of such a disk should reside (Novikov & Thorne 1973). Yet, we now 
know that when the disk dynamics is driven by magneto-rotational effects the source of free energy is not the gradient in the 
radial distribution of angular momentum but rather the gradient in angular velocity^ ; that is, the source of free-energy is the 
shear of the congruence of circular geodesies. Thus, MRI-mediated, turbulent angular momentum transport goes on unabated 
at the radius of marginal stability for cold geodesic orbits (Araya-Gochez 2002) and the flow within need not preserve specific 
angular momentum nor loose its relativistic Keplerian rotation profile in spite of a non-trivial radial velocity. 

We will therefore assume a Keplerian angular velocity profile throughout: r > rmb- This assumption is consistent with 
the models of PWF and Popham & Gammie (1999)-indeed, strong deviations from Keplerian rotation only occur as fluid 
elements approach the photon radius rph-and it is also consistent with global 3D simulations of non-radiative MHD accretion 
onto non-spinning black holes (Hawley & Krolik 2001, Hawlcy & Balbus 2002) in spite of partial pressure support for hot 
inner tori. 

Inside rms, the flow may not have time to cool significantly and advection of entropy will become progressively more 
important as the horizon, r+, is approached. Entirely advective, accretion flows (i.e., non-radiative, hydrodynamical flows) 
arc reckoned to posses a large, positive Bernoulli function B = q + J^dpp/g (Abramowicz et al. 1978, Blandford & Begelman 
2003) which makes the ultimate fate of accreting fluid particles a theme of controversy; e.g., ADAF (Narayan & Yi 1995) vs 
CDAF (Narayan et al. 2000) vs ADIOS (Blandford & Begelman 1999, 2003). The broad conclusion drawn from this ongoing 
debate is that fluids with large internal energies, e.g., strongly magnetized fluids, are not easily accreted onto gravitational 
wells. In the ultra-relativistic flows of interest to us, and for large spin parameter a, the value of the relativistic enthalpy, 
Q = (p + e+p), can bo non-trivial g/ p > 2 (with magnetic energy folded into the internal energy e), implying a high likelihood 
for large Bernoulli function as well. We adopt the view that for a fluid with such high enthalpy, the dynamical boundaries set 
by the circular orbits of cold (point) particles in the Kerr geometry are inadequate. 



3 MAGNETO-ROTATIONALLY-INDUCED RELATIVISTIC DISK DYNAMICS 

The MRI enables accretion through vigorous, turbulent transport of angular momentum. It also constitutes the key process to 
tap the free-energy available in a flow endowed with differential rotation. In the Newtonian picture, this is essentially a local, 
co-moving instability as the interplay of inertial accelerations with the elastic coupling of fluid elements creates an unstable 
situation to the redistribution of specific angular momentum, l(r), if the angular velocity profile decreases monotonically with 
radius. This criterion is in opposition to the (non-magnetic) Rayleigh criterion for stability of a differentially rotating fluid: 
r~^ dr£^ = «^ > 0, where k is the frequency of radial epicyclic motions. Indeed, without the elastic coupling provided by the 
bending of field lines, such inertial forces-namely, the shear(tide) and the coriolis terms-induce radial epicyclic motions while 
preserving specific angular momentum. 

Incompressible MRI initiated turbulence has a characteristic "parallel" scale for fastest growing modes corresponding to 
fe|| = k - 1b — O/iij^jf, and a normalized growth rate cr = A = ^dinrlnO, associated with the shear parameter, a.k.a. the Oort 
A "constant" . Although slower growing modes occur at larger scales, we take the fastest growing modes to be the dynamically 
predominant eddies given that the radial flow velocity is non-negligible in rmb < r < rms. Thus, the scale associated with 
the MRI is generally smaller than the vertical scale height of the disk by a factor 0[Vj^-^J cn] and no large-scale effects are 
expected. Indeed, in the weak field limit one may construct a dispersion relation quite independently of the large-scale field 
topology: highly sub-thermal fields, Vj^^f/cn <^ 1, guarantee that the instability is truly local and incompressible (see below). 
When the field is more moderate, fejj"^ will approach the disk's pressure scale height Ti ~ cn/^- and this imposes a lower 
bound on fc|| when the field is mostly vertical. On the other hand, the lowest bound on fcy in a toroidal field geometry is 1/r 
which would correspond to a supra-thermal field in the Newtonian picture (Fogfizzo & Tagger 1995, hereafter FT). Because 
of the possibility of large-scale, global disk dynamics, these modes are the focus of our analysis below. 



^ More generally, Balbus (2002) has found that for a magneto-rotatlonally-driven disk, the proper replacement of the classical H0iland 
stability criteria for rotational and convective motions involves the gradients in angular velocity and temperature instead of angular 
momentum and entropy. 
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3.1 Compressible Non-asymmetric MRI modes 

In the 2D regime of fastest growth, fee ^ fcr S> fe<^, non-axisymmetric "horizontal" displacement modes of a toroidal field have 
the following dispersion relation (Araya-Gochez 2002) 

- {(A + l)ql + x'} a' +Aql{ql + 4A} = 0, (1) 

where all frequencies are normalized to the rotation rate, A = idinrlnn is the Oort A "constant", = 4(1 + A) is the 
squared of the epicyclic frequency, and = (k • VAif)/f2 is a frequency related to the component of the wave vector along 
the field (in velocity units). 
A is defined through 

where F is the adiabatic index, and G = v'^^^/c^. This definition entrains an important anisotropy constraint on the Lagrangian 
displacement vector field ^ : v^^^ (kx ■ = — (k • ^), as well as the compressibility characteristics of non-axisymmetric 
modes (as a function of the component of Lagrangian displacement along the field, Foglizzo and Tagger 1995): 

^ = (l-A)(-ifc||^||) (3) 

where A denotes the Lagrangian (co-moving) perturbation as opposed to the Eulerian perturbation 5 (recall the uon-rclativistic 
relation A = 6 + ^ ■ 'V). Evidently, the compressibility of the modes is imprint on the deviations of A from unity. A word 
of caution here regards the peculiar behavior of the fluid when the internal stress has heat conduction characteristics, e.g. 
photon or neutrino pressure; sec §3.3. From Eq [3], one roads that the degree of compression of the modes gets stronger with 
the field strength and, naturally, with a softer equation of state. Note that setting A = 1 in the dispersion relation Eq [1] 
yields the incompressible, local variant of the MRI. 

Ftom the dispersion relation Eq [1], one finds the parallel wave- numbers of fastest growth to satisfy (see also Foglizzo 
1995) 

g| = -2A+(i^)x{-2AAi| -A(2 + A), where D ^ 1 + (l-A)A + ^1 + (1 - A)A, (4) 

while the expression in the curly brackets corresponds to the square of the growth rate. 



3.2 General Relativistic Effects in the Cowling Limit 

General relativity modifies this Newtonian results in that it introduces large-scale efi^ects when shear begins to overwhelm 
the coriolis terms, i.e., as Vph is approached (see Eqs [32] of Araya-Gochez 2002). Although the hole is likely to be born in an 

excited state, e.g., away from its asymptotic, Kerr geometry; we work in the Cowling limit, 5g = 0, and, to zeroth order, use 
the standard form of the Kerr metric in the equatorial plane (Boyer-Lindquist coordinates): 

ds^ = -^dt'^ + r^Aidf - ujdtf + ^dr'^, (5) 

with uj = 2a/ Ar"^ the rate of frame dragging by the hole and where the metric functions of the radial BLF coordinate are 
written as relativistic corrections (e.g. Novikov and Thorne 1973): 

.Asl-h — -1-2— , and X>sl---|-— , 

in normalized geometrical units (c = G = A/,^,^ = 1)- 

For circular, incompressible, geodesic flow, the dispersion relation near a rotating hole is given by (Araya-Gochez 2002) 

(7a)^- [24 +47^01 (C±-fP)] (7a)^ + g| [4 - 4 { §7^1) }] = (6) 

where 7 is the red-shift factor, © = and C± = 1 — | ± corresponds to the C function of Novikov & Thorne (1973) 
for pro-grade (-I-) orbits. Regarding the electromagnetic stress, the only difference with the non-relativistic, incompressible 
analog is that the Alfven speed is now weighted by the relativistic enthalpy of the fluid (Araya-Gochez 2002) 

gvl, ^ iF'-'F^.. (7) 

Factoring out the orbital frequency u''' /u* = Cl± = ±(r^/^ ± o)~^ and with 717 = Q± &, one has 

*^-[2C+xya^ + C[C+4A]=0 (8) 
where 

A = - { §7^0} and xl = 47' (C± - |o) 
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Figure 1. Normalized growth rate, (T, as a function of radius for several values of the spin parameter a (from Araya-Gochez 2002). 



denote the normalized shear parameter and co- moving epicycle frequency (i x corresponds to the epicycle frequency as 
measured at asymptotic infinity). Thus, with the proper generalizations of the epicycle frequency and shear parameter, the 
local dispersion relation is identical with the Newtonian case in the limit of no fluid compression, c.f. Eqs [1&3] with A ^ 1. 

Next, using the relation 7^ = (1 ± a/r^^^y C^^ for circular, geodesic flow (Novikov and Thorne 1973), one finds the 
fastest growing modes will conform with 



g. = 1 — 



J_ ^4 
16 X 



1± 



3 V 



(9) 



which remains finite and close to the Newtonian value — A(2 + A) = 15/16 for all radii outside of the ISCO. 

Recall that the marginally stable orbit (a.k.a. the ISCO), rms, corresponds to the root of x± = 0. The radius of the 
circular photon orbit, rph, is where C± = and the event horizon, r+, happens at the outer root of "P = 0. For any value of 
the rotation parameter a : rms > r-^h > ^ph > r+. Our statement on large-scale effects follows from these remarks and from 
direct inspection of Eq [9]: — > 0"*" as r — > r'^^, i.e. the most unstable MRI modes are go to large scales as the photon orbit 
is approached. Figs 1 & 2 show the general trends for the normalized growth rate and wave-numbers as functions of radius 
and spin parameter. 

Two notable caveats of the relativistic approach involve the global disk structure and the compressibility of the modes. 
We discuss these in turn. 

Curvature effects and radial field structure will indeed modify the dispersion relation^ (see, e.g.. Curry & Pudritz 1995, 
Ogilvie & Pringle 1996) but these effects are unlikely to modify the gross properties of the fastest growing modes unless 

radial stratification plays a destabilizing role. If such stratification were magnetically-driven; it would require an unlikely 
supra-thermal field v^^^^ « w*^ = f{Q. — lo), where f = rA/\/V is the radius of gyration for the physical velocity in the locally 
non-rotating frame (Bardeen, Press & Tcukolsky 1972). Furthermore, we note here and correct our previous comment (Araya- 
Gochez 2002): — > 0"'" is merely a statement on the length-scale of the fastest growing modes, fe|j~^ ~ 0[f\, regardless of field 
strength. 



^ the interested reader is invited to compare the relativistic, local version of the equations of motion for Eqs [32] of Araya-Gochez 2002, 
with the Newtonian, global version Eqs [2.6 &; 2.7] of Ogilvie & Pringle 1996 
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Figure 2. Normalized parallel wavenumber, ij^, as a function of radius (in gravitational radii) for several values of the spin parameter 
a. Diamonds indicate the location of the marginally stable orbit, x = 0) ^^'^ triangles, the location of the marginally bound orbit (from 
Araya-G6chez2002). 

The compressibility of the modes alters the threshold of shear parameter where ^ 0"'" , but docs not fundamentally 
affect the conclusions drawn from Eq [9] either. For roughly equal neutrino and radiation pressures (see below), one can use an 
effective adiabatic index (Araya-Gochez & Vishniac 2002) in Eq [4] to anticipate that the effects of compressibility on toroidal 
modes is to increment the threshold of shear parameter where 0^ from —A = 2 to —A ^ 1 + 2D /{I + A); c.f. Eq [4]. 

Nevertheless, since for geodesic flow A oc T'/C± and C± ^ @ rph, the increase in shear threshold in this setting is rather 
inconsequential. 

3.3 Heat Conduction and Clump Scale 

A precise assessment of radiative diffusion effects in the general relativistic regime is hampered by tfie breakdown of a key 
assumption made to simplify the "linear poking" on the Faraday field tensor (Araya-Gochez 2002) : use of the enthalpy weighted 
specific four magnetic field in Eq [7]. However, just as in radiation pressure dominated fluids, one expects a neutrino stress 
dominated gas-with neutrinos semi-contained by the nucleon component through scattering-to comprise a rather peculiar 
MHD fluid. Indeed, when a fluid transitions into a radiative pressure regime, compressive modes will loose pressure support 
in a given range of wavo-rmmbers (Agol & Krolik 1998). Thus, the magnetic field is truly frozen only to the co- moving volume 
associated with the baryon component and compressive perturbations need not behave adiabatically. 

The half optical depth to scattering of a standard radiation-pressure-dominated a-disk is related to the scale height, H, 
and rotation rate, Q,, by Tdisk = —c/{2aAQ'H), This relation is explicitly sensitive only to the local nature of the cooling but is 
implicitly subject to a suitable vertical gradient of heat deposition (see, e.g., Krolik 1999). Therefore, these scalings also hold for 
a neutrino cooled Q-disk through a replacement in the source of opacity by neutrino scattering in the non-advcctivc accretion 
regime: 0.1 Mq sec~^ < M < 1 M© sec^^. On the other hand, in the Newtonian regime q^ ~ 0[1] and for magnetic angular 

momentum transport such that a ~ v^it/c^ i the depth through the MRI eddies is Teddy ~ ''^au/'^h ^ ''"disk ~ c/(2Ay^c„) S> 1. 
Thus, assuming an isotropic random walk for neutrinos, the diffusion time through these eddies, tdifi ~ is similar to the 
time scale for fastest MRI modes to develop tmiu — ^AqJ^.^. 

Radiative heat conduction will isotropize the dynamically dominant modes since neutrinos will diffuse through the smallest 
distance associated with the eddies, i.e. mostly in the vertical direction. In our problem, this will diminish the growth rate of 



6 R. A. Araya-Gochez 



the instability (Blacs & Socrates 2001) but not by much. Indeed, if neutrino stress were predominant, the horizontal regime 
of fastest growth (e.g. §3.1) would be inaccessible; yet, such is never the case since the pressure contributions from radiation, 
pairs and neutrinos, all have the same temperature dependence a/6 with the relative contributions varying only by internal 
degrees of freedom times particle statistics factors (with only one helicity state for neutrino pairs): 2 x 1,4 x 7/8, and 6 x 
7/8 respectively. Thus, pi, is never greater than about ~ prad where prad = ll/12oT'' includes the pressure from pairs and 
photons. 

Let's quantify these arguments in order to gauge the length scale of the clumps. At a very fundamental level, clump 
formation is intimately connected to the effects of radiative heat conduction out of compressive perturbations. This can 
be understood by writing down the polarization properties of the fastest growing modes in the horizontal regime: = 
— a/a & l^el -C |^r|, and reading from Eq [3] that at the linear stage of the instability there exists a converging flow 
toward the Lagrangian displacement node of the modes (see Fig 2 of Foglizzo and Tagger 1995). In a fluid with entirely elEistic 
(adiabatic) properties, the pressure perturbation associated with such compression will act as a restoring force to de-compress 
the fluid in the non-linear stage. On the other hand, when the fluid has radiative heat conduction properties on the scale of the 
density perturbations, no such restoring force persists on time scales longer than the inverse of the growth rate so the clumps 
will survive for longer times. Turner et al. (2001, 2003) report that in a standard radiation-pressure-dominated disk, when 
Prad ^Pgas — Pb 1 the non-linear outcome of the MRI is a porous medium with drastic density contrasts. Under nearly constant 
total pressure and temperature, the non-linear regime shows that density enhancements anti-correlate with azimuthal field 
domains-just as expected from the linear theory-and that turbulent eddies live for about a dynamical time scale while mass 
clumps are destroyed through collisions or by running through localized regions of shear on a similar time scale. Notably, the 
non-linear density contrasts may be quite large < pmax/pmin > ^ 0(10]. We make no attempt here to estimate Ap/p from the 
lineaj theory but rather take the view that when the length scale of the instability is suitable large, Amri — 0[f] (see below), 
a fraction of C'[l/27r] of the mass in the annulus rmb < r < Tms will reside in one massive clump such that the gravitational 
wave emission associated with the in-fall of this clump will dominate the power spectrum. 

Let us emphasize that the isotropic size of the clumps is determined by the condition that, on the time scale of the 
instability, neutrinos will diffuse through the clumps via nucleon scattering and/or absorption followed by re-emission. For 
simplicity, we again assume local balance of dissipation from the shear flow with cooling from neutrino emission and a constant 
supersonic Mach number for the disk flow: = v"'' jcji {M,^ = 2.5 is consistent with the results of PWF in the neutrino 
cooling regime). Under this assumption, we need not worry about detailed neutrino transport and, relaxing the condition 
(jg ~ 1, the size of the most massive clumps will correspond to 

Zciump — fAit/^'w ~ \/aM~^r. 

Thus, in the relativistic regime the length-scale of the clumps will be smaller than the scale associated with the fastest growing 
modes by a factor of 



4 A SIMPLE PHYSICAL MODEL FOR GRAVITATIONAL WAVE EMISSION 

The key to our argument for clumpy, near-hole hyper-accretion is based on the observation that the length-scale of the fastest 
growing modes may be quite large at the innermost disk boundary rmb- Since the radial velocity is non- negligible, one cannot 
demand "tightly fitting" eigenmodes on the circumference at fixed radius. Loosely speaking, the WKB modes live on ingoing 
spiral trajectories (see §4.2); but wo arc only concerned with their overall gross properties. 

The clumps form very rapidly, tj^^ ~ fio", and very close to the horizon. Since these arc very compact, l/r ~ O[0.1], 
we considered them as point particles in the order of magnitude estimate of gravitational wave emission below. Computation 
of accurate waveforms for the in-fall of a single clump constitutes a daunting task because of the non-trivial multipole 
contributions to the radiation field in the strong gravity regime near rmb. The contributions from higher than quadrupole will 
broaden the power spectrum of emitted waves (Zerilli 1970, Davis et al. 1971, hereafter DRPP) and recent results by Lousto 
& Price (1997, hereafter LP97) show that particles released from very near the marginally bound orbit of a non-spiiming 
hole-or more precisely near the maximum of the Zerilli potential-will maximize the locally radiated energy. 

As a "back of the envelope" estimate of the quadrupole r.m.s. strain, \h\ = -^/i^ -|- /i ^ > associated with the in-fall of 
a single clump, we simply assume that the clump is formed above rmb and that it survives for one rotation period before 
plunging into the hole. The transverse-traceless projection of the metric perturbation in the radiation zone, hj^ , is estimated 
by replacing the second time derivative of the mass quadrupole moment, by the kinetic energy associated with 

non-spherical motion; i.e., clump rotation as seen at asymptotic infinity (Thorne 1987). Thus 



\h\ = 14^1 = ^ 



dtJi,it-r) 



■^n%,, GSM 



For a 15 Mq black hole rotating at a = .98, and accreting 1 Mq sec~^, the mass in the clumps is 5M ~ M/f2 ~ 1.8_4 Mq, 
where n(rmb) — .4 = 5400 Hz; corresponding to a linear frequency f2(rmb)/7'' = 1720 Hz. This gives ~ 3.2_25 for a source 
at 27 Mpc. To bring this into the LIGO band, \h\cH — I.-21 one would need to integrate over N ~ l+s cycles or over 1700 
sec! Thus, the r.m.s. strain associated with in-fall of a single clump onto the hole is too small to be of any astrophysical 
importance at present. 

But this issue is more subtle than it appears at first sight. 
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4.1 Quasi-normal Ringing Waveforms: Single Excitation Event 

Since the clumps form in a strong gravity regime, the process of gravitational wave emission from in-fall at such a short 
range must be cast as an excitation of the infinite number of increasingly damped quasi-normal modes of oscillation of the 
background geometry (Leaver 1985). The often-drawn analogy for such quasi- normal ringing (QNR) of the hole to the ringing 
of a bell is fundamentally flawed in at least one respect: the black hole is not excited by the smashing of the clump as it "hits" 
the horizon. The hole is rather excited when the metric perturbation associated with the clump is "felt" by the background 
metric. The excitation event therefore constitutes a smooth process whereby in-fall of a clump from ~ Vmh and through re- 
serves as a source in the Teukolsky (1973) equation for small perturbations to the Kerr geometry (with appropriate boundary 
conditions at r+ and Too)- This is an important distinction with a great deal of relevance to the problem at hand since we 
need to gauge the "driving" of QNR modes in terms of an effective coupling from clump in-fall. 

Clumpy black hole accretion from an exterior disk will excite preferentially quadrupolc, bar-like mode perturbations of 
the geometry with spheroidal harmonic indices^ (^, w) = (2,2). An approximate analytical expression for the QNR frequency 
of this mode is given by Echeverria (1989): 

0,22 = [1 - 0.63(1 - a)^/i«] X {l + ^(1 - a)«/^°} . (10) 
The corresponding quality factor, a.k.a. Q-value, of the (2,2) mode 

, , _ 1 3fJ[a)22] _ , ^-9/20 /^s 

Q22(a)= - -2(1 -a) , (11) 

Z ^[0^22 J 

is > 2tt for Ocrit > 0.92 This represents a minimum value of a to regard the hole as a decent bar-like-mode oscillator. Below 
ttcrit the mode is damped too strongly to allow for any significant build up of the strain signal/spectrum from resonant 
accretion (see below). 

In the radiation zone and time domain, undriven QNR waveforms constitute circularly polarized, damped oscillations 
(Leaver 1985, Echeverria 1989) 

h{t) = h+- ih^ = ^S22{(t>,e,a) e-'^"^^*-'^) (12) 

where d is the luminosity distance to the source, S22{(t>,0,a) is a spin weighted, normalized (J rff2|6'22|^ = 1), spheroidal 
harmonic of coordinates 4> Sz 6 [in the Schwarzschild limit 822(4', 0, a ~* 0) — > l22(</', 9)]. The strain length amplitude of 
the waveform. Ho (in units of the hole's mass M), and the phase shift, if, depend on the initial conditions of the metric 
perturbation. 

Ho gauges the efficiency of the emission process to posit a fraction of the total rest-mass energy of the system into 
outgoing gravitational wave energy. For strong excitation events such as binary black hole coalescence, this fraction may be as 
large as 3% (see, e.g. Flanagan & Huges 1998). For the in-fall of a clump of mass SM <^ M, this fraction scales as £{SM)^ /M 
(DRPP), with a corresponding time domain strain amplitude (x ^/e j \fM . The precise value of e for a rotating hole is 
unknown, but it hovers on a few per cent for the perturbations of a Schwarzschild geometry induced by radial in-fall from 
infinity'* (DRPP). For non-axisymmetric perturbations, the energy released from clump in-fall may be considerably larger 
(Fryer, Holtz & Hughes 2002) and for nearly resonant, driven oscillations from hyper-accretion, we demonstrate below that 
the total amount of energy deposited into outgoing gravitational waves is only bounded from above by a factor of a few times 
the rest-mass-energy of a single clump! 

4.2 Collective Effects: Resonant Driving of QNR modes 

To assess the free-fall stage of the clumps, we recall from §2 that the cusp in effective potential occurs somewhere above rmb 
(depending on angular velocity profile) and that a relativistically hot accretion flow may posses a semi-Keplerian angular 
velocity profile, 0+ = (r''^^ + a)^^, down to this region, (PWF, Popham & Gammic 1998, Hawlcy & Krolik 2001). Parameter- 
izing the accretion flow through its Mach number and the angular momentum transport through a Shakura & Sunyaev "a" 
parameter (as in §3.3, we adopt M,^ = 2.5 & a = 0.1), when oS^.9 the wavelength associated with the fastest growing MRI 
modes becomes comparable with 27rf at a radius r > Tmb; e.g., when ~ .126^ao.i M~ 2 5- At higher spin parameter values, 
MRI modes go to comparably large scales at radii closer to rms- Thus, the annulus rms > r > rmb constitutes a relativistic 
locus for large-scale magneto-rotationally induced fluid dynamics at moderately large values of black hole spin: a > 0.9. 

Since the radial velocity is non-negligible, MRI wave-modes must be interpreted as WKB modes that live on in-spiral 
trajectories. For the sake of simplicity in the discussion below, presume that fluid elements undergo one full rotation in the 
annulus rms > r > rmb before "plunging in" from rmb- The integrated phase along this path can be interpreted with a skewed 
effective wavelength which goes to increasingly larger length-scales as fluid particles approach rmb. If we demand a tight flt 

Quadrupolc perturbations, 1 = 2, are reckoned to contribute about 90% of tlic energy budget for outgoing GW's while 1 = 4 modes 
add about 10%. Higher multipole contributions are negligible, see, e.g. DRPP, LP97. 

* Strictly speaking, the energy released in this case includes not just the ringing phase but also a small contribution from the gravitational 
Bremsstrahlung component, see, e.g., LP97. 
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Figure 3. Shown are the real part of the QNR frequency, ojqnr, and the driving frequency, 0;^^ = 2^-1- (r j^b ) 1 for quadrupole oscillations 
as functions of the hole's spin parameter a. Plotted also is half of the imaginary part of the QNR frequency offset by +0.5 from 
zero for ease of comparison with the difference among the other two frequencies. Most notable is the fact that throughout the range 
a > 0.95, |t.iJQNH, — ^di\ ^ ^TgNR which implies resonant driving of the QNR modes from hyper-accretion. 



of one whole, integrated WKB wavelength on the in-spiral trajectory between Vms and rmb, the mid Lagrangian displacement 
node, ^(r) = 0, represents a converging point for the flow in the linear stage of the instability as seen by a co-moving geodesic 
observer. This node subsequently develops into the mass over-density that we identify with a clump in the non-adiabatic stage 
(i.e., when neutrinos have diffused out of the compressive perturbation, §3.3). Since the scale of the mode is larger downstream 
of the displacement node, most of the mass in the over-density will be advected in from fluid particles that rush upstream 
toward the node on the instability time scale; a smaller fraction of the clump's mass flows downstream toward the node, from 
smaller scales upstream. 

When the effective WKB scale on the instability is suitably large, a reasonably large fraction of the mass in the annulus will 
reside in a single clump by the time it reaches rmb. The subsequential "free-fall" from rmb (presuming this radius corresponds 
to the cusp in effective potential; e.g., Kozlowski et al. 1978, Abramowicz et al. 1978), occurs on a dynamical time scale: 
tdyn = n^^(rmb). Furthermore, since /-pole modes couple to the /-multipole of the dynamical frequency (e.g., Ex. 3.7 of 
Rybicky & Lightman 1979), the driving frequency for quadrupole oscillations corresponds to twice the relativistic Keplerian 
frequency at rmb : i^dr = 2r2+(rmb). Remarkably, the difference between the driving frequency and the resonant (2,2) QNR 
frequency resides unthm a factor of the order of the damping rate throughout the range .9 lia S.99 (Fig 3); whereas the Q value 
of the mode hovers on 5.6 < Q < 33. 

This picture of hyper-accretion would portray an ideal scenario to drive the QNR modes resonantly were the clumps to 
arrive at rmb steadily, with a small spread around time intervals of ©[nijl^], and for long enough to build the strain amplitude 
up to a saturation point-or about 2Q cycles (see below). It is noteworthy that relativistic magneto-rotational dynamics (§2) 
indicates that mass over-densities form rapidly enough to feed the QNR modes of the hole efficiently. Yet, nature will add 
randomness to the process. 

The added stochasticity in the arrival times of the clumps is related to the seed density perturbations in the flow. The 
accretion flow outside rms is turbulent on small scales, ^turb ^ r, and these density perturbations act as seeds for clump 
formation as these enter the relativistic annulus. This fact notwithstanding, in the picture of clump formation outlined above 
the probability distributions for arrival times and clump mass are strongly affected by large scale magneto-rotational dynamics 
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Figure 4. The solid line indicates the augmented amplitude of QNR oscillations driven near resonance by clumpy hyper-accretion from 
a neutrino-cooled disk as a function of the hole's spin parameter a. Modulo a small spread caused by the stochasticity in time arrival of 
the clumps to r^i-, (see discussion in the text), the enhanced amplitude lies within the curves for Q(a) and 2Q(a) throughout the range 
a > .9. A brisk energy deposition rate oc {Hsat / Hq)^ ensues. 



near rmb- If the flow is steady, large scale dynamics should minimize the impact of small scale perturbations upstream and 
induce a correlation between the two distributions. Below, we simply assume that the probability distribution for arrival times 
will favor resonant driving of the modes with a moderate spread around the resonant frequency of arrival times (see further 
discussion on this point below). Lastly, note that the arguments above do not depend on in-spiral trajectories wrapping around 
fully on the annulus; since the MRI time scale is suitable fast, a small delay to reach the free-fall stage is all that is needed 
for MRI dynamics to act. 

Beyond an initial transient period where the amplitude of the oscillation builds up to a saturation point, driven QNR 
oscillations become undamped oscillations at the driving frequency 



h{t) = 



Ho 



2 



d 



iWdrrQNR 



Si 



~~d 



S22(<j}, 0, a)e 



(13) 



where 5t — tan^^ (LJdrrQNR/[i.^dr ~ '^qnhJ) is a slowly varying function of time and the saturation amplitude, -ffsat, represents 
an effective strain length (in units of the system's total mass M) set by the condition that the energy input from clump in-fall 
(cx Hsat) be balanced by the damping losses (cx H'^^^). When the driving frequency is near resonance, ffsat is limited primarily 
by the natural^ damping of the QNR mode, FgNR = 9[lj22], and secondarily by the relative difference of the driving frequency 
to the resonant frequency for quadrupole ringing, lliqnr = 5R[tJ22]. Furthermore, since tJdr ~ i^qnr for a relatively broad range 
of a values (Fig 3), the resonant frequency may be taken to be ojqnr. 

The energy flux carried by driven QNR waveforms benefits greatly from the build up in the amplitude of the strain, 
-ffsat, when compared to the ringing from a single excitation event Hoic.i., Eq [12]). Indeed, when driven near resonance the 
augmented amplitude is proportional to twice the Q value of the mode: -ffaatUdf=t^QNr( = 2Q_ffo (c.f. Eq [11]). Although the 



^ Quasi-normal modes are strongly damped because of the B.C. of purely ingoing radiation at tj^ which saps most of the energy emitted 
by the metric perturbation unless the waves are strongly beamed from relativistic motion (Sterl Phiney, Priv. Comm.). In this sense, the 
damping of the mode corresponds to radiation reaction with most of the energy being advected in by the hole. 



10 R. A. Araya-Gochez 



modes are driven precisely at resonance only for a ~ .98, note from Fig 3 that Iojqnr — i^drl < Tqnr throughout the range 
.90SaS.99. This sets a lower bound on the augmented strain amplitude for this range of a: Hsat/Hoi^Q (recall that the 
resonant denominator only looks like a Lorentzian very near the resonance). Fig 4 shows the exact value of the augmented 
strain amplitude for a > .9. Furthermore, adopting a very modest £o = .02 for the in-fall of a single clump (in accordance 
with the DRPP and LP97 results; see also discussion in §4.1), the augmented ofEcicncy factor in GW energy deposition from 
driven oscillations e = (-ffsat/^?o)^£o may be larger than unity. Depending on the hole's spin, e ~ 1. — 10 for .92 < a < .98 
(c.f. Fig 4). 

Next, let us relate the energy emitted by the sequential in-fall of N clumps of mass SM through r+ to the amplitude 
of the waveform in the time domain. Assume N to be small enough such that N5M < Minitiai but large enough to enable a 
sensible average of the RMS strain amplitude over the corresponding number of wave cycles. The energy flux averaged over 
N cycles is (Isaacson 1968) 



dE 
dAdt 



G 



— {h{t)-h'{t)). 



(14) 



With h(t) as given by Eq [13] and assuming no spin evolution, d = 0, the total energy radiated in gravitational waves results 
from integrating Eq [14] over an encompassing surface (recall J dn]S'22]^ = l)i ami over a time interval corresponding to N 

cycles of LUdr'- 



Eg-w = 



c 

8G 



NLOdr\H^B.t\ 



(15) 



On the other hand, an approximate estimate of the total GW energy deposition which accounts for a slight mass increment 
of the black hole without spin evolution is as follows. First, write the energy emitted during one driven cycle as follows 
5E = [c^] (eSM) (5M/Minitiai, and note that the clump mass 5M ~ MCl':^_^{r■mb) oc M (constant o and M). For AM = I15M < 
A/initial, replace the factor in parenthesis by the geometric mean of SM over the hyper-accreting epoch, 5M, and perform the 
integral of 5M/M{t) by writing 5M = Mdt. This yields 



ASGw = [c']£<5Mln(^^) 

V Aiinitial / 

where e may be larger than unity for black hole spins in excess of 



(16) 



5 DETECTABILITY CRITERION 

Rather than attempt to describe the complex time evolution of the gravitational waveform throughout the hyper-accreting 
stage, we will simply assess the detectability of the waves with the matched filtering technique by utilizing the above energy 
estimate, Eq [16], while focusing on a relatively short time interval such that the added mass to the hole is a small fraction 

of the initial mass. 

For definitcncss, we adopt Minitiai = 15 M© (f ~ 2147 Hz) and A M = 1 M© corresponding to a second of hyper-accretion 
at M = 1 Mq sec~^. Under the simplified model for clump mass probability distribution as outlined above (quasi-linear 
analysis), a large fraction of the mass in the locus r-mh < r < r^s at any given time resides in a clump by the time fluid 
particles reach rmb- We estimate the mass in this clump to be 6M = Mf2^^(rmb) which is a factor of 27r smaller than the 
total mass in the annulus for fully wrapping, in-spiral trajectories, §4.2. 

We follow Flanagan & Huges (1998) to evaluate the likelihood of detection of GW from driven QNR modes by the LIGO 
interferometers. The criterion for detectability of the wave train through the matched flltering technique is a non-trivial S/N 
ratio 

([S/N]^)= /dln^^HllH (17) 

J ''noise J 



where /inoise(w) (= /inoise(/)/v27r) reflects the detector spectral density of strain noise appropriately averaged for random 
incident orientations with respect to the LIGO interferometers. The characteristic strain amplitude of the wave train is given 
by 

,2 , ._2(l + zfdE 

hche.rH= '^2^2 —[{1 + Z)0J], (18) 

where d is the luminosity distance to the source, z the redshift, and lj the observed angular frequency. Defined this way, /ichar 
is not trivially related to the amplitude of the strain in the time domain, i/gat, nor to the same in the frequency domain. 
It is rather a gauge of the time integrated power spectrum weighted by the "mono-chromaticity" of the wave train over the 
observation time. 

With (1 -I- z)lv —> uj, the GW energy spectrum from driven QNR modes in the radiation zone. 
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ensues from a Fourier transform of the strain signal in the time domain as given in Eq [13] (compare this expression with Eq 
[2.38] of Flanagan & Huges 1998, and note the conforming notation used throughout). 

Let us approximate all quantities oc M by their geometric mean and assume no spin evolution during a small epoch 
of hyper-accretion T = 2-iTN/uJdr (with logN in the range 3-3.5). With H{u)) = 2i!/'sate"'* sin[iV7rAc<;/wdr]/Aa; (note that 
d = 6t = constant as well) , the characteristic strain amplitude of driven QNR waves turns out to be 



/i-char(w) 



{I + zf Hs,t ■ H*,, 



Nnu) sin (a; — b) 



X — b 



(20) 



where x — b = Nit{u! — a;dr)/aidr. 

The strength of the characteristic strain signal at the driving frequency w = 2Q+ (rmb) follows by setting the expression 
in the square parenthesis to [Niv]'^ and by utilizing Eq [15] to replace |i?Bat|^ in favor of the energy released during the 

hyper-accreting epoch AEqw as given by Eq [16] 



,2 I 

'^char 



2{l + z)\ NAEgw 



Trd 



UJdr 



G 



27rc 



2(1 + g) 
7rd 



eSMT In I 



VMinitial/ ' 



(21) 



For GRB 030329, we plug in z ^ .1685, d = 810Mpc {Ho = 70, Qb = .3, Qv = .7), Mnitiai = 15 M©, T = 1 sec, and 
5M = MO,'^^{rinb) — 1.83-4 and choose a = .98 {Hsa.t/Ho = 22) to obtain /ichar — 8.4_23 at the observed linear frequency 
/ = 1490Hz. Utilizing similar parameters for GRB980425 at d = 27Mpc, turns out a characteristic signal /ichar — 2.I6-21 at 
/ ~ 1741Hz. 



6 DISCUSSION 

We have laid out a physical-if simple minded-model for resonant driving of the quasi-normal ringing (QNR) wave modes of 
the Kerr geometry. A nascent black hole hyper-accreting at rates M ~ 1 Mq sec^^ from a neutrino cooled disk is reckoned to 
oscillate near resonance of its {I, m = 2,2) quadrupole QNR frequency due to the in-fall of compact mass over-densities from 
the cusp in effective potential on a dynamical time scale. This model is based on large-scale, magneto-rotationally-induced 
fluid dynamics in the ultra-relativistic region of the flow bounded from below by the marginally bound orbit radius: rmb 
(assumed to coincide with the aforementioned cusp). The exact location of the cusp will change the dynamical time scale 
for in-fall but not by much when the spin parameter a > .9. Furhermore, since the MRI time scale is suitable fast, a small 
delay to reach the free-fall stage is all that is needed for MRI dynamics to act. Heat conduction from neutrino diffusion out 
of mass over-densities on the time scale of the instability warrants that compressive perturbations behave non-adiabatically. 
This leads to sufficiently long-lived clumps in a highly intermittent medium. When the spin parameter of the hole exceeds 
90%, large-scale magneto-rotational wave modes will segregate compact matter over-densities from large-scale magnetic field 
domains. The subsequent in-fall of these clumps from rmb will drive the QNR modes of the geometry in resonant fashion if 
the clumps arrive to rmb steadily, with a small spread around time intervals of 0[f2^^], and for long enough to build the strain 
amplitude up to a saturation point. 

It is difficult to address the probability distribution for clump arrival times from the linear theory alone but one anticipates 
that it will involve a co-moving time scale of ©[fi^^] and the integrated, spin dependent WKB wavelength on the inspire 
trajectory toward rmb- We have argued that for large mass clumps, 5M ~ Mf2^^, the arrival times to rmb will not be entirely 
stochastic since "the bucket needs time to fill-up" , nor entirely coherent since the seed density perturbations are stochastic 
but large-scale MRI dynamics will work to minimize the impact of this initial small-scale randomness. If the arrival times were 
entirely stochastic, the enhancement in energy deposition would be roughly proportional to the Q value of the mode instead 
of to Q'^ and the corresponding signal would be smaller by about an order of magnitude. A detailed probability distribution of 
arrival times requires an accurate account of the radial velocity profile and of the effective, spin-dependent WKB wavelengths 
on in-spiral trajectories which is beyond the scope of this study. 

If the QNR modes arc fed resonantly for a few seconds of hypcr-accretion, the enhanced amplitude of the oscillations 
yields a very high rate of energy deposition into gravitational waves. Indeed, the integrated energy deposition is large enough 
to "evaporate" the equivalent of a factor of a few times the total rest- mass-energy of a single clump into gravitational waves. 
Yet, the characteristic frequency of emission resides on a less than optimal range with respect to the advanced LIGO noise 
curve: / ~ 8000 — 1600Hz for M ~ 3 — 15 Mq. On the other hand, the high quality and coherence of the expected signal 
should make this model a prime candidate for narrow-band LIGO searches. For gammar-ray bursts at low redshifts, this may 
ultimately lead to a signal strong enough to be seen by LIGO II. 

A 15 M0 black hole accreting at 1 M© sec~^ may come about from the death throes of a very massive star. If such an 
explosion is associated with a long GRB, and the broad properties of the prompt gamma-ray light curve, to inner engine 
activity; hyper-accretion onto this heavy weight, fast spinning hole should be associated with a second broad hump in the 
light curve. Wc anticipate that gravitational waves with the predicted properties may be present at this stage of the GRB. 
The proposed mechanism is thus a direct probe of the elusive inner engine of gamma-ray bursts. 
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